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ABSTRACT 

We test an analytic model for the two-point correlations of galaxy clusters in redshift 
space using the Hubble Volume N-body simulations. The correlation function of clus- 
ters shows no enhancement along the line of sight, due to the lack of any virialised 
structures in the cluster distribution. However, the distortion of the clustering pattern 
due to coherent bulk motions is clearly visible. The distribution of cluster peculiar 
motions is well described by a Gaussian, except in the extreme high velocity tails. 
The simulations produce a small but significant number of clusters with large peculiar 
motions. The form of the redshift space power spectrum is strongly influenced by er- 
rors in measured cluster redshifts in extant surveys. When these errors are taken into 
account, the model reproduces the power spectrum recovered from the simulation to 
an accuracy of 15% or better over a decade in wavenumber. We compare our analytic 
predictions with the power spectrum measured from the APM cluster redshift survey. 
The cluster power spectrum constrains the amplitude of density fluctuations, as mea- 
sured by the linear rms variance in spheres of radius 8/i^^Mpc, denoted by erg. When 
combined with the constraints on erg and the density parameter derived from the lo- 
cal abundance of clusters, we find a best fitting cold dark matter model with erg « 1.25 
and 17 w 0.2, for a power spectrum shape that matches that measured for galaxies. 
However, for the best fitting value of and given the value of Hubble's constant from 
recent measurements, the assumed shape of the power spectrum is incompatible with 
the most readily motivated predictions from the cold dark matter paradigm. 
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1 INTRODUCTION 

Rich clusters of galaxies are unique tracers of the large scale 
structure of the Universe. This is due to a number of rea- 
sons. First, clusters are the most massive virialised systems 
in place at any given epoch, occupying a special place at the 
head of the structure hierarchy. Second, although rich clus- 
ters are rare objects they are also bright, containing many 
luminous galaxies and emitting copious amounts of X-rays. 
Third, a cosmologically interesting volume of the universe 
can be surveyed much more rapidly using clusters than with 
galaxies, as a result of the huge difference in the number 
of redshifts that are required to be taken. Finally, the key 
reason for the utility of clusters is that it is far easier to in- 
terpret the observed properties of the cluster distribution in 
the context of theoretical models than is the case for galax- 
ies. 

There are three main statistical properties of the cluster 
distribution that have been used to place constraints upon 



the parameters of structure formation models. The first of 
these is the local abundance of rich clusters, which, as a con- 
sequence of their rarity, is extremely sensitive to the ampli- 
tude of density fluctuations on a scale that encloses the mass 
of a typical rich cluster before collapse (White, Efstathiou 
& Frenk 1993). The rate at which the abundance of clusters 
evolves with redshift further constrains the mass density of 
the universe (Eke, Cole, Frenk & Henry 1998; Blanchard 
et al. 1998). The second test is the distribution of cluster 
peculiar motions, which depends on the value of the mass 
density parameter SI, for models with similar mass fluctu- 
ation amplitudes (Croft & Efstathiou 1994b; Bahcall, Cen 
& Gramann 1994). Lcistly, the spatial correlations of clus- 
ters constrain both the shape and amplitude of the power 
spectrum of the underlying dark matter (Mo, Jing & White 
1996). 

The latter of these, the clustering of clusters, is by far 
the best studied and, at the same time, the most controver- 
sial. The first measurements of the spatial two-point correla- 
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tion function of Abell clusters demonstrated a much stronger 
clustering amplitude than that found for galaxies (Bahcall & 
Soneira 1983; Klypin & Kopylov 1983). Moreover, the clus- 
tering amplitude was found to increase significantly as the 
cluster number density decreased (Bahcall & West 1992). 
However, the exact correlation amplitude of clusters remains 
the subject of intense debate (Efstathiou et al. 1992; Miller 
et al. 1999) . The early redshift surveys drawn from the Abell 
catalogue showed a significant enhancement of the clustering 
signal along the line of sight (Bahcall & Soneira 1983; Post- 
man, Huchra & Geller 1992). Bahcall, Soneira & Burgett 
(1986) found that the clustering in the line of sight direction 
is consistent with a "velocity broadening" of 2000kms~^, 
which they interpreted as arising from a combination of pe- 
culiar motions and geometrical distortions of superclusters. 

Confronted by these results, several authors have ar- 
gued that the Abell catalogue is afflicted by the superposi- 
tion of clusters and that the clustering signal along the line of 
sight is artificial (Sutherland 1988; Sutherland & Efstathiou 
1991; see also Lucey 1983). This prompted the construc- 
tion of more objectively defined cluster catalogues drawn 
from machine-scanned survey plates with better calibrated 
photometry (APM: Dalton et al. 1992, 1994, 1997; Cosmos: 
Lumsden et al. 1992). The typical radius used to define clus- 
ters in the machine based catalogues is significantly smaller 
than that used by Abell, reducing the enhancement of clus- 
ter richness by projection effects. The clustering signal found 
in these more recent cluster redshift surveys does not display 
large enhancements along the line of sight; furthermore, the 
trend of increasing correlation amplitude with decreasing 
space density of clusters is weaker than that found for Abell 
clusters (Croft et al. 1997). A similar dependence of corre- 
lation length on cluster space density is found for the X-ray 
Bright Abell Cluster Sample, which is much less susceptible 
to line of sight projection effects than the optically selected 
Abell catalogue (Abadi, Lambas & Muriel 1998, Borgani 
et al. 1999). This result was recently confirmed using the 
REFLEX survey in a work by CoUins et al. (2000), who 
found that there is no significant dependence of the clus- 
tering amplitud on the limiting flux, or equivalently on the 
space density of galaxy clusters. 

Miller et al. (1999) counter these objections by point- 
ing out that the early redshift surveys of Abell clusters con- 
tained large fractions of low richness clusters (Abell rich- 
ness class R = Q), that were not intended to form complete 
samples for use in statistical analyses. Furthermore, many 
cluster positions were determined by a single galaxy red- 
shift. Miller et al. (1999) present the clustering analysis of 
a new redshift survey of Abell clusters with richness -R > 1, 
and with the majority of cluster positions determined us- 
ing several galaxy redshifts. The clustering signal along the 
line of sight is greatly reduced in the new redshift surveys 
compared with the Bahcall & Soneira (1983) results, and 
is comparable to the amount of distortion of the clustering 
pattern found for APM clusters (see Fig 5 of Miller et al. 
1999). The anisotropy is further reduced after the orienta- 
tion of two superclusters that are elongated along the line of 
sight is changed. Peacock & West (1992) also found that re- 
stricting attention to higher richness Abell clusters removed 
the strong radial anisotropy seen in the clustering measured 
in the earlier surveys. 

It is clearly important to establish exactly what influ- 



ence peculiar motions have on the inferred spatial distri- 
bution of rich clusters. In most previous theoretical stud- 
ies, redshift space distortions have either been ignored or 
treated in a very approximate fashion. Quite often, redshift 
space distortions are modelled by simply assuming a boost 
to the clustering signal measured in real space, as predicted 
by Kaiser (1987). This effect arises from coherent flows on 
large scales where linear perturbation theory is applicable. 
However, if the peculiar motions of clusters are significant, 
then a damping of the clustering signal is expected on small 
scales. The transition between these two extreme types of 
behaviour needs to be modelled. 

Computer simulations of structure formation through 
the gravitational amplification of small primordial fiuctua- 
tions have been used extensively to model the spatial dis- 
tribution of clusters (e.g. White et al. 1987; Bahcall & Cen 
1992; Croft & Efstathiou 1994a; Watanabe, Matsubara & 
Suto 1994; Eke et al. 1996a). These early studies do not 
reach a consensus on the predicted clustering of clusters in 
cold dark matter cosmologies. Part of the reason for this dis- 
crepancy is due differences in the way in which clusters are 
identified in the simulations (Eke et al. 1996a). A further 
issue is the relatively small simulation volumes used and the 
small numbers of clusters analysed. 

Recently, it has become possible to simulate much larger 
volumes than were used in these earlier studies, with suffi- 
cient resolution to allow the reliable extraction of massive 
dark matter haloes that can be identified as rich clusters 
(Governato et al. 1999; Colberg et al. 2000). In this paper, we 
analyse the redshift space clustering of massive dark matter 
haloes in the Hubble Volume simulations, the largest cosmo- 
logi cal simulations to date, which are described in Section 



2.1 



This extends the comparison carried out by Colberg et 
al. (2000), who measured clustering in real space and com- 
pared the results with the predictions of analytic models. 
We outline the analytic model that we empl oy t o predict 



clustering in redshift space in Sections 2.2 and 2.4 



The dis- 
tribution of cluster peculiar velocities in the simulations, an 
important ingredient of the analytic model for redshift space 
clustering, is analysed in Section 2.3. The predictions of the 
analytic model are confronted with measurements of clus- 
tering in the APM Cluster redshift survey in Section ^ We 
compare to cluster power spectrum data directly rather than 
to estimates of the correlation length. This avoids uncertain- 
ties introduced by the method used to derive a correlation 
length from the data. Moreover, the cluster power spectrum 
on large scales has the same shape as the power spectrum of 
the dark matter, as demonstrated in real space by Colberg 
et al. (2000), and, as we show in Section |2.4| , is also the case 
under certain conditions in redshift space. We discuss our 
results and the constraints on cosmological parameters from 
the power spectrum of clusters in Section ^ , 



2 THEORETICAL PREDICTIONS FOR THE 
CLUSTER POWER SPECTRUM 

2.1 N-body simulations 

We compare analytic predictions of the statistical proper- 
ties of cluster samples with measurements made from the 
Virgo Consortium's "Hubble Volume" simulations. The sim- 
ulations follow the evolution of Cold Dark Matter density 



© 0000 RAS, MNRAS 000, 000-000 



Cluster correlations in redshift space 3 



fluctuations in two cosmologies: rCDM (with cosmologi- 
cal parameters = 1, a power spectrum shape defined 
by r = 0.21, following the parameterisation given by Efs- 
tathiou, Bond & White (1992), and a rms linear variance 
on a scale of 8/i"^Mpc of as = 0.6) and ACDM (with 
Qo = 0.3, a cosmological constant Aqc^ / {3Hq) = 0.7, a 
power spectrum described by an effective shape parameter of 
r = 0.17 and erg = 0.9). The huge volume of the simulations 
(8/i"^Gpc^ for rCDM and 27/i"^Gpc^ for ACDM) and the 
large number of particles employed (10^) allow cluster statis- 
tics to be studied with unprecedented accuracy (Colberg et 
al. 2000; Jenkins et al. 2001). Dark matter haloes are iden- 
tified using a friends-of-friends algorithm with a standard 
linking length (see Jenkins et al. 2001). The halo peculiar 
velocity is the peculiar motion of the centre of mass. 

2.2 Clustering in real space 

In this section we review the formalism employed to model 
the power spectrum of clusters using positions measured in 
real space, i.e. ignoring any distortion to the power spectrum 
arising from the peculiar motions of clusters. For further 
details, we refer the reader to the more complete discussions 
of this framework given, for example, by Mo & White (1996), 
Mo, Jing & White (1996), Borgani et al. (1997), Colberg 
et al (2000) and Moscardini et al (2000). Moscardini et al. 
(2000) also consider the evolution of clustering along the 
observer's past light cone, an effect that we shall ignore for 
the relatively shallow observational sample studied in this 
paper. 

Throughout this paper, we consider cluster samples that 
are defined by a characteristic spatial separation, dc, or 
equivalently, by a space density n, where dc = Such 
a sample is constructed by first ranking the clusters in order 
of mass, and then, starting from the most massive cluster, in- 
cluding progressively less massive clusters until the required 
space density is achieved. 

We assume that on large scales, the real space power 
spectrum of the cluster sample, Pc{k), can be related to the 
power spectrum of the underlying dark matter distribution, 
P{k), by an effective bias factor, 6eff , that is independent of 
scale: 

P.(fc) = blsP{k). (1) 

Colberg et al (2000) use the Hubble Volume N-body sim- 
ulations to demonstrate that this is an excellent approxi- 
mation over a decade in wavenumber, 0.01/iMpc~^ < k < 
O.l/iMpc^^. The task of computing the real space power 
spectrum of clusters can therefore be broken down into two 
steps: (i) The calculation of the appropriate power spectrum 
of the dark matter distribution, P{k), taking into account 
non-linear evolution of density fluctuations, (ii) The compu- 
tation of the linear bias factor, 6cff, for clusters of a given 
abundance. 

The first stage in the calculation is carried out using the 
prescription for transforming a linear theory power spectrum 
into a non-linear power spectrum described by Peacock & 
Dodds (1996). The transformation depends upon the cos- 
mological parameters Q and A, and upon the epoch or nor- 
malisation of the linear theory power spectrum. The formula 
given by Peacock & Dodds agrees well with the non-linear 
evolution found in N-body simulations (Jenkins et al. 1998) . 



The effective bias is computed by taking a weighted 
average of the bias, b{M), for haloes of mass M over the 
cluster sample under consideration 

TOO ^(^)dn(M)^^ 
TOO dn(M),,. ' 



where Afum is the lower mass limit that defines the sample 
and dn/dM is the space density of halos in the mass inter- 
val M to M + 5M (e.g. Mo & White 1996; Governato et 
al. 1999). We adopt the analytic form for the mass function 
of dark matter haloes proposed by Sheth, Mo & Tormen 
(2001 - hereafter SMT). These authors put forward a mod- 
ification to the theory of Press & Schechter (1974) in which 
the collapse of dark matter haloes is followed using an ellip- 
soidal rather than spherical model. The SMT mass function 
agrees well with the results of N-body simulations, although 
the most significant improvements over Press-Schechter the- 
ory are realised for lower mass haloes than we consider in 
this paper (SMT; Jenkins et al. 2001). Following the the- 
ory developed by Mo & White (1996), SMT also derive an 
expression for the bias factor of dark matter haloes (their 
equation 8) which we adopt in our calculations. 

As reported by Colberg et al. (2000), the SMT formu- 
lae for the mass function and halo bias factor predict an 
effective bias that is in good agreement with the results ob- 
tained from the Hubble Volume simulations. The analytic 
predictions for the real space power spectrum of clusters 
with dc — 30.9h~^Mpc are shown by the solid lines in the 
upper panels of Fig. ^ The discrepancy is largest for the 
ACDM model, in which case the analytic prediction for the 
real space power spectrum is 12% higher than the measure- 
ment from the simulation. 

2.3 Cluster peculiar velocities 

The gravitationally induced peculiar motions of clusters dis- 
tort the pattern of clustering if cluster redshifts are used to 
infer their spatial distribution. The distribution of cluster 
peculiar motions is therefore a key ingredient in the theo- 
retical prediction of clustering in redshift space, as discussed 
in the next section. 

In Fig. 0(a) and (b), the histograms show the distribu- 
tion of line of sight pairwise peculiar velocities, f{vi2), mea- 
sured for clusters in the Hubble Volume simulations with 
dc — 30.9/i~^Mpc; (a) shows the distribution for the rCDM 
simulation and (b) shows the results for ACDM. The smooth 
curves show various analytic distributions plotted with the 
pairwise velocity dispersion that is measured in the simula- 
tions {\/ < > = 532kms~^ for clusters in the rCDM sim- 
ulation and v^< > = 434kms"^ for ACDM); the solid 
line shows a Gaussian distribution, the dashed line shows an 
exponential distribution in |fi2|'^'^^, the dotted line shows 
an exponential distribution in |wi2|'^''* and the dot-dashed 
line shows an exponential distribution in |tii2| (full details of 
the functional forms of these distributions may be found in 
Padilla & Baugh 2001). The bulk of the distribution of clus- 
ter peculiar velocities is adequately described by a Gaussian. 
However, Fig. |l| illustrates that this is not the case for the 
high velocity tails of the distributions. Moreover, the best fit- 
ting distribution is not the same in each simulation. Padilla 
& Baugh (2001) show that the shape of the high velocity 
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Figure 1. The histograms show the distribution of line of sight pairwise pecuUar velocities of clusters with dc = 30.9h~^Mpc in the 
Hubble Volume simulations: (a) shows the results for tCDM and (b) for ACDM, as indicated in the legend. The lower panels show the 
distributions when a single cluster rms redshift error of 5v = 500kms~^ is included. The smooth curves show theoretical distributions 
plotted with the same variance in pairwise peculiar velocity that is measured from the simulations. In (a) and (b), solid lines show 
a Gaussian distribution, dashed lines show an exponential distribution in |?;i2p''^, dotted lines are for an exponential distribution in 
|^'12|''^*, and the dot-dashed lines show an exponential distribution in \vi2\. In the lower panels, the variance includes the redshift error 
described above, added in quadrature to the variance measured in the simulation. In (c) and (d) only the Gaussian distribution is plotted. 



tail depends upon the degree of non-linear evolution of the 
density fluctuations. The distribution of peculiar velocities 
is therefore of interest in itself, being sensitive to the cos- 
mological parameters Q and erg, and to the power spectrum 
of density fluctuations (Croft & Efstathiou 1994b). These 
issues are explored in more detail using the Hubble Vol- 
ume simulations and other N-body simulations by Padilla 
& Baugh (2001). 

The histograms in the lower panels of Fig. ^ show the 
distribution of peculiar velocities after incorporating a single 
cluster rms redshift measurement error of Sv ~ 500kms~^, 
appropriate for the Abell and APM cluster redshift sur- 
veys (Efstathiou et al. 1992). The resulting distributions are 
much broader and closer to Gaussian, as shown by the cor- 



responding solid lines; in these cases the variance is given by 
the sum in quadrature of the variance measured in the simu- 
lation and the rms redshift error stated above. The redshift 
error dominates over the variance expected for gravitation- 
ally induced peculiar velocities for the models that we con- 
sider in this paper. We therefore make the approximation 
in subsequent calculations when comparing to APM data 
that the variance in the pairw ise peculiar velocity, including 
errors, is fixed at a/ < 1^12 > ~ 850kms~^, regardless of the 
model being studied. This in t urn imp lies a single particle 
rms peculiar velocity of ct = < ^12 >/V2 ~ 600kms~^. 

2.4 Clustering in redshift space - power spectrum 
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2.4-1 Analytic model 



2-4-2 Comparison with simulation results 



The distortion of the power spectrum measured in redshift 
space generally displays two forms. On large scales, the am- 
plitude of the power spectrum is boosted due to coherent 
inflows into overdense regions and outflows from underdense 
volumes (Kaiser 1987). On small scales, randomised motions 
within virialised structures cause a damping of power (e.g. 
Peacock 1999). Simple models have been developed that de- 
scribe the transition between this large and small scale be- 
haviour (Peacock & Dodds 1994, Cole, Fisher & Weinberg 
1995). These schemes have been shown to work reasonably 
well for dark matter in N-body simulations (Cole, Fisher & 
Weinberg 1995; Hoyle et al. 1999). In this section we test 
whether such models provide an accurate description of the 
redshift space power spectrum of clusters of galaxies; this is 
necessary as we do not expect to flnd virialised structures 
in the cluster distribution. 

For scales that are still evolving according to linear per- 
turbation theory, the power spectrum in redshift space is 
given by 

2\2 



P<f(fc,M)=Pc(fc) 1+/3M 



(3) 



where Pc is the cluster power spectrum in real space, as 
defined by equation ^, Pc{k) is the cluster power spectrum 
in redshift space, /3 = /(n)/6(,fi (/(fi) is the logarithmic 
derivative of the fiuctuation growth rate) and n is the cosine 
of the angle between the wavevector k and the line of sight 
(Kaiser 1987; see also the discussion of this result in Cole, 
Fisher & Weinberg 1994). 

Heuristic schemes have been put forward that extend 
this model for the redshift space power spectrum down to 
small scales to include the effects of a random velocity dis- 
persion, under the assumption that the velocities are uncor- 
related with the density field: 



P<f(fc,Ai) = P,{k) (l -H/3m') D(k^Ia^)- 



(4) 



For the case of a Gaussian distributed velocity dispersion, 

(5) 



DikjiOv) ~ exp(— fc^/i^(T^/2) 



whilst for an exponential distribution (Ballinger, Peacock & 
Heavens 1996), 



D(kfiav) ~ 



1 



(6) 



1 + {kiiG^Y /2 

The spherically averaged form of equation ^ for a Gaussian 
velocity dispersion is given by Peacock & Dodds (1994): 

P:,{k)^G{l3,y)P^i{k), (7) 

where the function G{I3, y) is given by 

G{P,y) =^^[3/?^+4/?y^+4y4] 



CXp(-I/ ) \r,2 



[f3^(3 + 2y^)+4l3y% 



(8) 



where /3 = /(r2)/6eff and y = fc(T„/(100kms~^Mpc"^). Er- 
rors in the determination of the cluster redshifts can be in- 
corporated into this model by adding the redshift error in 
quadrature to the rms peculiar velocity to redefine o"„. On 
the scales that we consider, there is effectively no difference 
in the distortion to the power spectrum when a Gaussian 
or exponential distribution of velocity dispersion is adopted 
(see Fig. 3 of Ballinger, Heavens & Peacock 1996). 



The analytic prediction for the spherically averaged cluster 
power spectrum in redshift space is compared to the mea- 
surements from the Hubble Volume simulations in Fig. ^ 
In the simulations, the clustering pattern in redshift space 
is obtained by displacing clusters along the x-axis by an 
amount 5x = v^/ Ho, where is the x-component of a clus- 
ter's peculiar motion. The power spectrum of clusters in the 
Hubble Volume simulation is computed using the technique 
described in detail by Jenkins et al. (1998), and is equiva- 
lent to employing a very high resolution grid to perform the 
Fast Fourier Transform (FFT). Therefore, the scheme used 
to assign clusters to the FFT grid does not influence the 
recovered power spectrum. 

The lower panels in Fig. |^ show the ratio of the red- 
shift space to real space power spectrum as a function of 
wavenumber. The dashed lines show the outcome of divid- 
ing the analytic prediction for the real space spectrum into 
the prediction for the redshift space power spectrum, calcu- 
lated including the effects of cluster redshift errors. At small 
wavenumbers (large scales) the ratios are in excellent agree- 
ment with the expectations from equation ^, which, using 
the approximation /(fi) ~ ^^'^ predicts a ratio of 1.22 for 
clusters in rCDM and 1.15 in the ACDM simulation. Note 
the latter value changes by less than 1% if the weak influence 
of a non-zero cosmological constant on /(f2) is taken into 
account (Lahav et al. 1991). At high wavenumbers (small 
scales) we find that there is no damping of the power spec- 
trum measured in redshift space when redshift errors are 
ignored. The analytic model reproduces the boost in power 
on large scales but predicts too much damping in the power 
on small scales. If an error is included in the cluster redshifts, 
then the model predicts the same form of distortion found 
in the simulations. The upper panels of Fig. ^ show that 
in this case, the analytic model for the redshift space power 
spectrum agrees with the simulation results to 15% or better 
over the wavenumber range —2 < log(fe//iMpc~'^) < —1. 

Further insight into the form of the distortion of the red- 
shift space power spectrum caused by cluster motions can 
be obtained by plotting the power spectrum as a function of 
wavenumber perpendicular (fcj_) and parallel (fc||) to the line 

of sight, P{k±, fc||), where k = \Jk\ ~\- and ^ = k^^/k. We 

plot P(A;_L,fc||) for dc = 30.9/i"^Mpc clusters in the rCDM 
Hubble Volume simulation in Fig. pi The smooth lines show 
analytic predictions. The light solid lines show the real space 
power spectrum. In the upper panel, the heavy solid lines 
show the redshift space power spectrum. In the lower panel 
the heavy lines show the power spectrum in redshift space 
including the effects of redshift errors in the determination of 
cluster positions. Two contour levels are plotted; the inner- 
most sets of contours show P{k±,k\^) — 4.9 x 10'*/i~^Mpc^ 
and the outermost set show P{k±, fc||) = 1.2 x lO^/i'^Mpc^. 
On large scales, the analytic predictions are in excellent 
agreement with the measurements from the simulations, in 
real space and in redshift space. In redshift space, the power 
is enhanced on large scales in the fc|| direction, displacing 
the contour of fixed power to higher wavenumbers. On in- 
termediate scales, the distortion of the power spectrum is 
extremely sensitive to how well cluster redshifts are mea- 
sured. The magnitude of the error estimated in the redshifts 
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Figure 2. The upper two panels show the real and redshift space power spectra for clusters with dc = 30.9ft.~^Mpc; the left panels 
show the results for tCDM and the right panels show results for ACDM. The lines show analytic predictions and the symbols show 
measurements made directly from the simulations. Real space quantities are shown by solid lines and squares. The redshift space power 
spectrum measured in the simulations is shown by the circles. The dashed lines and triangles show the redshift space power spectrum 
after incorporating a redshift error into the determination of cluster positions. The lower panels show the ratio between redshift space 
and real space power spectra as a function of wavenumber. Again, points show the ratio for the measurements in the simulation. The 
circles show the ratio when only peculiar motions are considered, the triangles show the ratio when redshift errors are also included. 
The dashed line shows the analytic prediction of the ratio, when cluster redshift errors are included and should be compared with the 
triangles. The error bars are set by the number of modes per bin in wavenumber in the simulation. 



of APM and Abell survey clusters dominates over the distor- 
tion due to the peculiar motions of the clusters. Therefore, 
on these scales, the boost in power given by equation ^ is a 
poor description of the power spectrum. Qualitatively, the 
analytic model reproduces the form of the redshift space dis- 
tortion measured in the simulation. However, due to a small 
discrepancy on intermediate scales between the predicted 
real space power spectrum and the measurement from the 
simulation (of around 10% in the amplitude of P(fc)), the 
contours do not coincide on this plot. 

Finally, we compare the shapes of the redshift space 
power spectrum of clusters and dark matter in the rCDM 



Hubble Volume simulation in Fig ^. The results for clusters 
include errors in the cluster redshift determination, as dis- 
cussed above. This plot illustrates that the redshift space 
power spectrum predicted by the model is in very good 
agreement with that measured for the dark matter. The ef- 
fective bias measured in redshift space, as deduced by taking 
the ratio of cluster and mass power spectra in redshift space, 
is somewhat lower than the bias measured in real space. 
However, the bias measured in redshift space is still inde- 
pendent of scale. A measurement of the cluster power spec- 
trum on large scales in redshift space would therefore yield 
the shape of the mass power spectrum in redshift space. 
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Figure 3. The power spectrum as a function of wavenum- 
ber parallel A:|| and perpendicular k± to the line of sight. The 
smooth lines show theoretical predictions, the other lines show 
measurements from the rCDM Hubble Volume simulation. The 
light lines in both panels show the power spectrum in real space. 
In the upper panel, the heavy lines show the power spectrum 
when the peculiar motions of clusters are taken into account. In 
the lower panel, redshift errors in the cluster positions are in- 
cluded. The inner most set of contours show the wavenumbers for 
which the power is 4.9 X 10''/i~^Mpc'^; the outermost set show 
P(fcx,fc||) = 1.2 X 10*/j-^Mpc^ 

2.5 Clustering in redshift space - correlation 
function 

A key statistic for cluster samples is the two point correla- 
tion function ^, measured as a function of cluster separation 
parallel to the line of sight, tt, and perpendicular to the line 
of sight, a. Efstathiou et al. (1992) used this statistic to 
argue that cluster samples drawn from the Abell catalogue 
are contaminated by projection effects, leading to a spuri- 
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Figure 4. A comparison of real and redshift space power spectra 
for dc = 30.9/i~^Mpc clusters and for dark matter in the rCDM 
simulation. Solid symbols and lines show real space quantities and 
open symbols and dashed lines show redshift space quantities. 
Squares show power spectra for the dark matter and circles show 
cluster power spectra. Redshift errors have been included in the 
redshift space power spectrum of clusters. The solid lines in the 
upper panel show the fits to the real-space CDM model for the 
mass and clusters, and the dashed lines show the same fits for 
redshift space. The lower panel shows the effective bias obtained 
by taking the square root of the ratio between cluster and mass 
power spectra, both in real space (solid points and lines) and in 
redshift space (open points, dashed line). 

ous enhancement of the clustering signal along the line of 
sight. Previously, the two-point correlation function ^(cr, vr) 
has been studied for more abundant clusters, using much 
smaller volume simulations than are considered in this pa- 
per (e.g. Eke et al. 1996a). The Hubble Volume simulations 
can be used to resolve once and for all the issue of exactly 
how much anisotropy is expected in the two point correla- 
tion function from redshift space distortions alone. 

In Fig ^ we show ^(ct, tt) for clusters defined by dc = 
30.9/i~^Mpc in the Hubble Volume simulations. The left 
hand panels show correlation functions measured in the 
rCDM simulation and the right hand panels show results 
for ACDM. The upper row of Fig. ^ shows the correlation 
function measured in real space and the lower row shows the 
redshift space correlation function. Note that in this Fig., 
we do not include any redshift errors. The contour levels are 
given in the figure caption; the thick contour shows ^ = 1. 
As expected, the contours do not show any distortion when 
the correlation function is measured in real space. However, 
when the effects of peculiar motions are included, we find an 
apparent enhancement in the clustering signal perpendicular 
to the line of sight. This is a result of the flattening of the 
contours in the it direction due to coherent flows in the clus- 
ter distribution. There is no evidence for any enhancement 
of the clustering signal along the line of sight. As a test, we 
have also measured ^((t, tt) for the dark matter in the Hub- 
ble Volume simulations, and in this case we do flnd a strong 
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Figure 5. The two-point correlation function tt) plotted in bins parallel (tt) and perpendicular (cr) to the line of sight for dc = 
30.9/i^^Mpc clusters in the tCDM (left-hand panels) and ACDM (right-hand panels) Hubble Volume simulations. The upper row 
shows the correlation function in real space and the lower row shows the redshift space correlation function. Contour levels are at 
^ = 3, 2, 1, 0.8, 0.6, 0.4, 0.2. The thick contour shows $ = 1. 



enhancement of the clustering amplitude in the tt direction 
on small scales. The relatively large separation of particles 
on the initial grid in the simulations (3/i~^Mpc in the case 
of ACDM) is therefore not an issue. The explanation of our 
result for clusters is that virialised structures have not had 
time to form in the cluster distribution. 

The correlation function measured in the rCDM simu- 
lation is shown on an expanded scale in Fig. ^ In the lower 
panel, we include a rms cluster redshift error of 500kms~^, 
in addition to the peculiar motions. When cluster redshift 
errors are included, a clear boost is evident in the amplitude 
of clustering parallel to the line of sight. It is interesting to 
note that the flattening of the contours due to coherent flows 
is no longer apparent; the magnitude of the cluster redshift 
errors is sufficient to obscure this effect. It can also be seen 
that the clustering amplitud along the a direction shortens 
noticeably. This can be explained as follows: the effect of 



the redshift errors on the correlation function at fixed val- 
ues of (T, will be that of a smoothing function, which will take 
power from the scales corresponding to the higher correla- 
tions and transfer it to those of lower correlation. Therefore 
the value of ^((J, tt = 0) will be smaller after the redshift er- 
rors are included, and the correlation lenght in the direction 
perpendicular to the line of sight will be shorter. 



3 COMPARISON OF THE THEORETICAL 
MODEL WITH DATA 

3.1 Clustering data 

In this paper we compare the predictions of CDM mod- 
els with the power spectrum measured from the APM clus- 
ter redshift survey. The data we consider are the variance 
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Figure 6. The two-point correlation function measured in red- 
shift space for clusters in the rCDM simulation. The same contour 
levels used in Fig. ^ are plotted, with the thick contours showing 
^ = 1. In panel (a), gravitationally induced peculiar motions are 
considered. In (b), a rms cluster redshift error of 500kms~^ is 
also included. 



of counts in cells measured by Gaztanaga, Croft & Dalton 
(1995) and the power spectrum measured by Tadros, Ef- 
stathiou & Dalton (1998). Both measurements were made 
using sample B of the APM cluster redshift survey, as de- 
fined by Dalton et al. (1994). 

Before confronting the model predictions with these 
data, we first compare the two measurements with one an- 
other. The variance of counts in cells is given by an integral 
over the power spectrum multiplied by the Fourier trans- 
form of the window used to smooth the density field (e.g. 
Peacock 1999) 



27r2 



dkk'^P{k)W^{kR). 



For a spherical top hat smoothing window of radius R, 
3 



W{kR) 



(kRy 



■ (sin(fci?) - kRcos{kR)) . 
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Figure 7. Testing the recovery of the power spectrum from 
the variance. The solid line shows a linear theory CDM power 
spectrum with erg = 1 and F = 0.5. The points show the power 
spectra recovered by applying equations |ll] and ^ to the variance 
computed from the linear theory power spectrum {as given by 
equation The different symbols denote the results obtained 
for different assumptions about the slope of the power spectrum. 



The integral is reasonably sharply peaked around a charac- 
teristic wavenumber for a particular smoothing scale. There- 
fore, we make the approximation that the variance measured 
in a sphere of radius R can be related to the power spectrum 
at a specified wavenumber (see Peacock 1991 who gives the 
result for a Gaussian smoothing window): 

a^{R)=A\k,B), (11) 

where A'^{k) = fc='P(fc)/27r^ If we consider power law spec- 
tra, P{k) oc fc", then the effective wavenumber is defined 
as: 



fccff = ^ [9/(n)l^/("+3) ^ 



(12) 



Here n is the logarithmic slope of the power spectrum at 
fceff, and the function I{n) is defined by (see also Gaztanaga 
1995): 



/(n) = 



dxx" 



(sin X — X cos x) 



(13) 



We demonstrate the accuracy of this transformation in Fig. 
^ The points show the power spectra recovered by applying 
equations ^ and ^ to the variance computed by inte- 
grating over the power spectrum shown by the solid line, 
which is our input spectrum (using equation The differ- 
ent symbol types delineate the results obtained when dif- 
ferent assumptions are made about the shape of the power 
spectrum. The most accurate answers are obtained when 
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Figure 8. (a) Estimates of the power spectrum of APM clus- 
ters inferred by applying equations [l| and to the variance of 
counts in cells measured by Gaztaiiaga et al. (1995). The differ- 
ent symbols show the results for different assumptions about the 
slope of the power spectrum: the solid triangles use the slope of 
the variance data directly, the upwards pointing triangles are for 
a fixed value of n = — 1 and the downwards pointing triangles 
are for n = —0.5. (b) A comparison of the power spectrum in- 
ferred using the counts-in-cells measurements of Gaztanaga et al. 
(1995) (filled triangles) with the direct measurement of the power 
spectrum of APM clusters by Tadros et al. (1998) (circles). 



the shape of the power spectrum is inferred directly from 
the local slope of the variance. In this case, the recovered 
spectrum is at most 10% below the correct value. The loca- 
tion of the turnover is reproduced when fixed values of the 
spectral index n axe used; however, the amplitude of the re- 
covered spectrum can differ from that of the true spectrum 
by a factor of two, on scales larger than the turnover. 

In Fig ^(a), we show estimates of the power spectrum 
made by applying equations and |l2| to the variance of 
counts in cells for APM clusters taken from Figure 3 of 
Gaztaiiaga et al. (1995). The different symbols show the re- 
sults for different assumptions about the logarithmic slope of 
the power spectrum: the filled triangles show the power spec- 
trum when the slope n is estimated directly from the mea- 
sured variance; the upwards pointing, open triangles show 
the results for a fixed value of n = — 1 and the downwards 
pointing triangles show the case where n — —0.5. The power 
spectrum recovered in this way from the variance is robust 
to reasonable variations in the value of the spectral index 
adopted in equation |l^. 

We compare the power spectrum estimated from the 
variance of counts in cells, using the value of n inferred from 
the data, with the power spectrum measured directly in red- 



shift space by Tadros et al. (1998) in Fig ^h). Over most 
of the range considered, the Tadros et al. (1998) result is 
approximately 50% higher in amplitude than that obtained 
by Gaztaiiaga et al. (1995) Furthermore, though both mea- 
surements agree within the errors on the largest scales, the 
suggestion of a turnover seen in the Tadros et al. (1998) 
power spectrum is not as apparent in the power spectrum 
inferred from the variance data. In both papers the same 
sample of clusters was considered, with the same cuts ap- 
plied in redshift. Both approaches make assumptions that 
could be partly responsible for the discrepancy. Tadros et 
al. (1998) need to specify a radial weighting scheme to ob- 
tain a minimum variance estimate of the power spectrum. 
The weighting scheme depends upon the power spectrum, 
but in practice the recovered spectrum is fairly insensitive 
to the value chosen. In addition, a random distribution of 
points with the same selection function as the cluster data 
is used to model the effects of the survey geometry on the 
recovered power spectrum. The measured power spectrum 
is sensitive to the level of smoothing applied to the selection 
function, though the tests presented by Tadros et al. (1998) 
suggest that the turnover found in the power spectrum is 
robust to changes in the smoothing. The Gaztanaga et al. 
(1995) approach does not require an estimate of the survey 
selection function (see Efstathiou et al. 1990), but does as- 
sume that the fluctuations are Gaussian. However, relaxing 
this assumption does not change the results significantly. 
Furthermore, the estimators used for the power spectrum 
and for the variance of counts in cells could respond in dif- 
ferent ways to uncertainties in the mean density of clusters 
and to the irregular boundary of the APM survey. 



3.2 Model parameters 

We now compare theoretical predictions for the power spec- 
trum of clusters in redshift space with measurements of the 
power spectrum of the APM cluster redshift survey. We 
make a number of assumptions in our model: 

(i) We assume that primordial density fluctuations are 
Gaussian. The bulk of the available evidence suggests that 
any non-Gaussianity found in the distribution of large scale 
structure today is consistent with the gravitational evolu- 
tion of an initially Gaussian density field (Moore et al. 1992; 
Gaztaiiaga 1994; Canavezes et al. 1998; Hoyle, Szapudi & 
Baugh 2000; Szapudi et al. 2000). Relaxing this assump- 
tion would affect the model predictions for the abundance 
of clusters and therefore, through changing the bias factor, 
the clustering of clusters (Robinson, Gawiser & Silk 2000). 

(ii) The dark matter in the universe is assumed to be cold 
dark matter. 

(iii) We assume that the shape of the dark matter power 
spectrum is the same as that measured for the galaxy 
power spectrum on the scales that we consider. Local bias- 
ing schemes yield an asymptotically constant bias on large 
scales between fluctuations in the galaxy and mass distribu- 
tions (Coles 1993; see Cole et al. 1998 for realisations of this 
process). Observational support for this assumption is de- 
rived from the comparison of the power spectra of different 
types of object. Peacock & Dodds (1994) found that on large 
scales, the power spectra of optical galaxies, radio galaxies 
and clusters are consistent with a single underlying power 
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spectrum after applying different, constant relative bias fac- 
tors to the individual measurements. Tadros et al. (1998) 
also demonstrate that the power spectra of APM clusters 
and galaxies have the same shape on large scales and are 
related by a constant relative bias factor. The connection 
with the shape of the dark matter spectrum is made theo- 
retically. Colberg et al. (2000) show that in real space the 
power spectrum of clusters is the same shape as the power 
spectrum of the underlying dark matter. In Section 



2.4 



have demonstrated that the same conclusions hold in red- 
shift space, when cluster redshift errors are included. 

The shape of the real space power spectrum of galaxies is 
well determined (Maddox et al. 1990; Baugh & Efstathiou 
1993, 1994; Peacock & Dodds 1996; Gaztanaga & Baugh 
1998). We adopt values of F = 0.2, 0.25, consistent with the 
recent determinations by Efstathiou & Moody (2001) and 
Eisenstein & Zaldarriaga (2001). 

(iv) We apply the methodology outlined in Section ^ and 
assume that the construction of a cluster sample specified by 
a space density is equivalent to taking all clusters above some 
mass threshold. To compare with sample B drawn from the 
APM cluster survey, we adopt a value of dc = 30.9/i~^Mpc. 

Specifically, our model has two parameters: 

(i) The cosmological density parameter, Qq- We consider 
values in the range ilo = 0.01 — 1, with the condition that the 
model universe is spatially flat, i.e. if < 1, then we adopt a 
cosmological constant. A, such that n-|-Ac^/(3-ff^) = 1. This 
is in agreement with the location of the first Doppler peak 
reported by the BoomeranG and Maxima cosmic microwave 
background experiments (Balbi et al. 2000; de Bernardis et 
al. 2000). 

(ii) The amplitude of density fiuctuations, as specified by 
the linear rms variance in spheres of radius 8ft~^Mpc, ug. 

The theoretical models are computed for a mean cluster red- 
shift of z = 0.083, which is appropriate for APM sample B. 

When comparing the model predictions with the data, 
we use the APM power spectrum results obtained for a cos- 
mology with f2 = 1. Tadros et al. (1998) also compute the 
power spectrum of APM clusters in a background cosmol- 
ogy defined by Qo = 0.2, Aoc^/{3Hq) = 0.8, and recover 
the same shape of power spectrum, but with an amplitude 
that is approximately 25% higher. Therefore we make a sys- 
tematic error when comparing data derived assuming Q — 1 
with a model in which a different value of il is adopted. To 
partially compensate for this, we take a pessimistic view of 
the errors on our theoretical predictions, adopting the level 
of the discrepancy found between the model and the ACDM 
simulation results of 15%. 



3.3 Results 

The theoretical models are assessed by computing a figure 
of merit: 



X 



= (Pmodcl(fc) - Pdata(fc))' / (aL^k) + aLdcl (fc) ) (14) 



where the sum is over all the data points plotted in Fig 
^b) and we take crniodGi(fc) = 0.15Pmodei(fc) (see Robinson 
2000). We ignore any covariance between measurements of 
the power spectrum at different wavenumbers. 



The results are shown in Fig. ^ The shaded regions 
show models that are within Ax^ < 1,4, 9 of the best fit- 
ting model when compared to the cluster power spectrum 
data indicated in the legend on each panel. The legend also 
gives the value of the power spectrum shape parameter, F, 
used in the models. In Fig. ^, the left hand panels show the 
results of the comparison to the Tadros et al. (1998) power 
spectrum data, and the right hand panels show the results 
using the power spectrum inferred from the Gaztanaga et al. 
(1995) measurement of the variance of APM clusters. The 
top row of the Fig. is for models with F — 0.2, the bottom 
row is for F = 0.25. We also plot the same Ax^ contours 
for the constraint on the value of as derived from the abun- 
dance of hot X-ray clusters in the local universe by Eke, 
Cole & Frenk (1996b) (shown by the solid, dashed and dot- 
ted lines). The error on erg quoted by these authors is 8%. 
The Eke et al. results agree with those from a recent reanal- 
ysis of the current cluster X-ray temperature data and of 
the theoretical framework of this constraint carried out by 
Pierpaoli, Scott & White (2000). The almost horizontal lines 
show the models in which clusters with a mean separation of 
dc = 30.9/i~^Mpc have effective bias factors of bcs — 1.5,2 
and 3. 

The Ax^ contours derived from the cluster power spec- 
trum constraint are almost parallel to the Q axis, particu- 
larly for n > 0.2. Thus the amplitude of the cluster power 
spectrum in the models is determined largely by the value 
of erg (see Mo, Jing & White 1996). The best fitting mod- 
els with F — 0.25 have a lower P^r degree of freedom 
than the F = 0.2 models. The difference in the amplitude of 
the Tadros et al. (1998) and Gaztanaga et al. (1995) data 
is readily apparent from the shift in the location of the 
Ax^ contours. The models that come closest to matching 
the Gaztanaga et al. (1995) data tend to have an effective 
cluster bias factor of befi ~ 2, whilst the best fitting models 
for the Tadros et al. (1998) data tend to have bes ~ 1.5. 

The constraints on model parameters derived from the 
power spectrum data are much broader than those obtained 
from the abundance of rich clusters. Moreover, the shapes 
of the two sets of contours are different. A rough fit to 
the middle of the Ax^ < 1 contour for the Gaztanaga 
et al. (1995) data, using models with F = 0.25, yields 
erg — j^~'' i*+'' 0'^(i~'^t))^ whereas the cluster abundance con- 
straint scaling is approximately as cc Q^"'^^. 

We illustrate the consequences of an error in the space 
density of APM clusters in Fig. |l^. The estimated space 
density is somewhat sensitive to the way in which the clus- 
ter selection function is normalised (Efstathiou et al. 1992). 
Furthermore, any projection effects that persist in the ma- 
chine constructed catalogues could be responsible for mov- 
ing poorer clusters into the sample. We therefore plot in 
Fig. |l^ the constraints on model parameters assuming that 
dc — 34/i~^Mpc; this represents a 10% error in dc, corre- 
sponding to the space density of APM clusters being 30% 
lower than assumed in Fig. ^. When the space density of the 
sample decreases, the minimum mass threshold that defines 
the sample increases, leading to larger effective bias parame- 
ters. The Ax^ contours therefore shift down to lower values 
of erg. 

Finally, we consider the constraints on the model pa- 
rameters that are obtained when the two datasets, the clus- 
ter power spectrum and local abundance of hot X-ray clus- 
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Figure 9. The constraints on the model parameters o-g and f2o from APM cluster power spectrum data. The shaded regions show 
contours of = li 4 and 9. The almost horizontal lines show parameter combinations for which the effective cluster bias is fegjf = 1.5, 2 
and 3. The left hand panels show the outcome of comparing the models to the Tadros et al. (1998) power spectrum data, the right hand 
panels are obtained using the power spectrum inferred from the Gaztaiiaga et al. (1995) variance of counts in cells. The upper row is for 
models with V = 0.2, the lower row for F = 0.25. The lines show Ax^ contours for the constraint on erg and Ho from the local abundance 
of hot X-ray clusters and are reproduced in each panel (Eke, Cole & Frenk 1996b). 



ters, are combined. In Fig. ^ we plot the Ax^ contours 
after adding the values from comparing the model to the 
power spectrum data and to the observed cluster abundance. 
This operation assumes that the two datasets are indepen- 
dent. Due to the smaller errors, the cluster abundance con- 
straint has the largest influence on the resulting contours. 
The vertical line on each panel shows the value of Q. that is 
consistent with the chosen F, given the recent determina- 



tion of Bubble's constant as Hq = 72 ± 8kms~^Mpc~^ by 
Freedman et al. (2001), and assuming that F = Qh, where 
Ho ~ 100/ikms~^Mpc~^. The dotted lines indicate the range 
of values of allowed following this prescription when the 
1(7 errors on Bubble's constant are taken into account. The 
combined dataset favours low values of n, with ag = 1-1.5 
for the Gaztahaga et al. measurement and ag = 1.6-2.3 for 
the Tadros et al. power spectrum. In both cases, the shape of 
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Figure 10. The constraints on model parameters when the 
abundance of APM clusters is reduced by 30%, so that dc = 
34/t~^Mpc. The lines and shading are the same as for Fig ^ 

the power spectrum expected in the most easily motivated 
cold dark matter model, given the best fitting value of fl is 
somewhat discrepant with the shape that we have assumed 
for the power spectrum; this disagreement is significant for 
the Tadros et al. data. 



4 DISCUSSION 

The main goal of this paper was to establish the degree of 
distortion expected in the clustering pattern of rich clus- 
ters when viewed in redshift space. For this purpose, we 
undertook high precision measurements of the clustering of 
massive dark matter haloes, with the same space density 
as observational samples of clusters, using the largest ex- 
tant cosmological simulations, the Virgo Consortium's Hub- 
ble Volume simulations. 

The rms peculiar motions of clusters in the Hubble Vol- 
ume simulations are modest, of the order 300-350kms~^ for 
a single cluster. Examples of pairwise velocity dispersions 
of the order of several thousand kilometers per second are 
found in the simulations, but are relatively rare. We find no 
evidence for any enhancement of the correlation amplitude 
along the line of sight arising from peculiar motions. The 
main type of distortion to the clustering pattern in redshift 
space is a fiattening of the contours of ^(cr, tt) due to coher- 
ent motions of clusters. The shape of the ^{a, tt) = 1 contour 
is remarkably similar to that found by Miller et al. (1999) 
for a recent survey of f? > 1 Abell clusters (see their Fig. 
7; we note however that lower amplitude contours in their 
Fig. still show some enhancement along the line of sight). If 
a redshift error is assigned to each cluster in the simulation 
(with a single cluster rms of 500kms~^), then we reproduce 
the form of the distortion seen in measurements of ^(""1"") 
from the APM cluster redshift survey, confirming that this 
survey is largely free from projection effects. 
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Figure 11. The best fitting model parameters when the con- 
straints from the cluster power spectrum and abundance measure- 
ments are combined. The contours show A(Xp^j.j ~^^abun^ = li 4 
and 9. The top panel shows the Ax^ contours when the Tadros et 
al. (1995) P{k) is used and the bottom panel is for the Gaztafiaga 
et al. (1995) data. A model power spectrum shape of F = 0.25 is 
adopted and the model clusters are defined by dc = 30.9/i~^Mpc. 
The vertical line shows the value of Q that is consistent with the 
adopted value of F, the CDM definition of F = Qh and the mea- 
surement of Hubble's constant by Freedman et al. (2001). The 
vertical dotted lines show the range of Q permitted within the la 
errors quoted on Hubble's constant. 



Similar conclusions were reached in a study of the 
redshift space distortions in the clustering of galaxies and 
groups in the Updated Zwicky Catalogue by Padilla et al. 
(2001). In this case, the groups are identified in three di- 
mensions rather than in projection. A different set of is- 
sues need to be addressed when tuning a three dimensional 
group-finding algorithm, but superposition of groups along 
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the line of sight should only occur for systems with high 
internal velocity dispersions. Padilla et al. (2001) found a 
strong enhancement in the line of sight clustering for galax- 
ies. No such feature was detected in the correlation func- 
tion of groups. These results were compared with predictions 
from N-body simulations, and the same form of distortion 
was found. 

The analytic model outlined in Section ^provides an ac- 
curate description of the redshift space power spectrum of 
clusters measured in Hubble Volume simulations, when er- 
rors in the cluster redshifts are included. We then compared 
the model predictions to the power spectrum measured from 
the APM cluster redshift survey in order to constrain the 
cosmological parameters Q, and cts • This is a more direct ap- 
proach than a comparison with the cluster correlation length 
and thus avoids uncertainties regarding the way in which the 
correlation length is derived from the measured correlation 
function. 

We have chosen to focus our attention on the APM 
cluster redshift survey, as this is the largest volume survey 
available that has been constructed from survey plates us- 
ing a well specified, automated procedure. Miller & Batuski 
(2001) measure the power spectrum of a sample of i? > 1 
Abell clusters that does not display large enhancements in 
the line of sight clustering and which covers a larger volume 
than the APM survey. These authors find no evidence for 
a turnover in the cluster power spectrum and probe larger 
scales than the measurements of Gaztanaga et al. (1995) and 
Tadros et al. (1998). However, the space density of clusters 
in Miller & Batuski's sample changes by a factor of two in 
a northern extension of the survey, which contributes much 
of the signal to the power spectrum measurement on large 
scales. Therefore more detailed modelling of the observa- 
tional sample is required than we have attempted in this 
paper in order to make a realistic comparison with these 
data. Furthermore, the integrity of the Abell catalogue re- 
mains open to question, even for R>1 clusters. Van Haar- 
lem, Frenk & White (1997) demonstrated that mock Abell 
cluster catalogues constructed from N-body simulations suf- 
fer from significant incompleteness for 7? > 1 clusters. Power 
spectrum measurements have also been made using X-ray se- 
lected samples of clusters (e.g. Zandivarez, Abadi & Lambas 
2001, Schuecker et al. 2001). Such samples have the appeal 
of greatly reducing projection effects as the X-ray emission is 
dominated by the high density core of a cluster's dark mat- 
ter halo. However, flux limited X-ray surveys mix clusters of 
different richness, so more careful modelling of the cluster 
selection is required to make robust theoretical predictions 
(Borgani et al. 1999; Moscardini et al. 2000). 

We assume that the dark matter power spectrum has 
the shape measured for the galaxy power spectrum on large 
scales. This is a reasonable approximation if galaxy forma- 
tion is a local process (Coles 1993; Cole et al. 1998). The 
model parameters that we vary to flt the cluster power spec- 
trum data are then the cosmological density parameter 
and the fluctuation amplitude as- Within this scheme, the 
cluster power spectrum does not provide any constraint on 
Q, (see Mo, Jing & White 1996). The constraints on erg are 
fairly broad as a result of the relatively large errors on the 
measured cluster power spectrum. This situation can be im- 
proved if additional information is used. Tadros et al. (1998) 
compare the power spectrum of galaxies and clusters in the 



APM survey and find that they have the same shape and 
can be related by a relative bias of b ~ 1.8. If one assumes 
that APM galaxies are essentially unbiased tracers of the 
dark matter on large scales (Gaztanaga 1994), then this ad- 
ditional constraint can be used to exclude models in which 
APM-like clusters have an effective bias that is very differ- 
ent from b ~ 1.8. However, we caution against the use of 
measurements of clustering in redshift space to restrict the 
value of an effective bias parameter in real space (see Fig. 
0)- 

The range of acceptable model parameters is tightened 
considerably if the cluster power spectrum data is combined 
with measurements of the local abundance of rich clusters. 
The cluster abundance measurements themselves constrain 
the parameter combination aaOP'^ (White, Efstathiou & 
Frenk 1993). The constraint on these parameters from the 
power spectrum data has a different dependence on f2, thus 
allowing the degeneracy to be lifted. The best model param- 
eters obtained using the combined data sets are as « 1.25 
and ~ 0.2. The best fltting value of Q. is inconsistent with 
the value suggested by our choice of F, if we use the most 
readily motivated prescription for setting the shape of the 
power spectrum in cold dark matter models, F = ^Ih, and 
adopt a recent measurement of Hubble's constant by Freed- 
man et al. (2001). The discrepancy is severest for the Tadros 
et al. (1995) power spectrum data (see the discussion of this 
data in Gawiser & Silk 1998). Whilst values of the shape 
parameter F < Q.h can be motivated physically by postulat- 
ing a large baryon fraction or decaying neutrinos (Eisenstein 
& Hu 1998 ; White, Gelmini & Silk 1995), it is difficult to 
see how a value for F that is much larger than could be 
accommodated. However, this situation could be alleviated 
if the abundance of APM clusters has been overestimated 
and a larger value of dc turns out to be more appropriate. 

The cluster power spectrum data will improve signifi- 
cantly in the near future upon completion of the 2dF and 
SDSS galaxy redshift surveys (CoUess 1999; York et al. 
2000). Group catalogues will be constructed in three dimen- 
sions and the large volumes covered by the surveys mean 
that these catalogues will contain large numbers of clusters. 
Moreover, a comparison between the properties of clusters 
selected in projection and in redshift space will permit fine 
tuning of the two dimensional algorithms, that can then be 
applied to the deeper parent photometric catalogues. 
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